2 results
A playground for compressible natural convection with a nearly uniform density
- Thierry Alboussière, Jezabel Curbelo, Fabien Dubuffet, Stéphane Labrosse, Yanick Ricard
-
- Journal:
- Journal of Fluid Mechanics / Volume 940 / 10 June 2022
- Published online by Cambridge University Press:
- 05 April 2022, A9
-
- Article
- Export citation
-
In the quest to understand the basic universal features of compressible convection, one would like to disentangle genuine consequences of compression from spatial variations of transport properties. For instance, one may choose to consider a fluid with uniform dynamic viscosity, but, then, compressible effects will generate a density gradient and consequently the kinematic viscosity will not be uniform. In the present work, we consider a very peculiar equation of state, whereby entropy is solely dependent on density, so that a nearly isentropic fluid domain is nearly isochoric. Within this class of equations of state, there is a thermal adiabatic gradient and a key property of compressible convection is still present, namely its capacity to viscously dissipate a large fraction of the thermal energy involved, of the order of the well-named dissipation number. In the anelastic approximation, under the assumption of an infinite Prandtl number, the number of governing parameters can be brought down to two, the Rayleigh number and the dissipation number. This framework is proposed as a playground for compressible convection, an opportunity to extend the vast corpus of theoretical analyses on the Oberbeck–Boussinesq equations regarding stability, bifurcations or the determination of upper bounds for the turbulent heat transfer. Here, in a two-dimensional geometry, we concentrate on the structure of numerical solutions. For all Rayleigh numbers, a change in the vertical temperature profile is observed in the range of dissipation number between $0$ and less than $0.4$, associated with the weakening of ascending plumes. For larger dissipation numbers, the heat flux dependence on this number is found to be well predicted by Malkus's model of critical layers. For dissipation numbers of order unity, and large Rayleigh numbers, dissipation becomes related to the entropy heat flux at each depth, so that the vertical dissipation profile can be predicted, and so does the total ratio of dissipation to convective heat flux.
Numerical solutions of compressible convection with an infinite Prandtl number: comparison of the anelastic and anelastic liquid models with the exact equations
- Jezabel Curbelo, Lucia Duarte, Thierry Alboussière, Fabien Dubuffet, Stéphane Labrosse, Yanick Ricard
-
- Journal:
- Journal of Fluid Mechanics / Volume 873 / 25 August 2019
- Published online by Cambridge University Press:
- 26 June 2019, pp. 646-687
-
- Article
- Export citation
-
We developed a numerical method for the set of equations governing fully compressible convection in the limit of infinite Prandtl numbers. Reduced models have also been analysed, such as the anelastic approximation and the anelastic liquid approximation. The tests of our numerical schemes against self-consistent criteria have shown that our numerical simulations are consistent from the point of view of energy dissipation, heat transfer and entropy budget. The equation of state of an ideal gas has been considered in this work. Specific effects arising because of the compressibility of the fluid are studied, like the scaling of viscous dissipation and the scaling of the heat flux contribution due to the mechanical power exerted by viscous forces. We analysed the solutions obtained with each model (fully compressible model, anelastic and anelastic liquid approximations) in a wide range of dimensionless parameters and determined the errors induced by each approximation with respect to the fully compressible solutions. Based on a rationale on the development of the thermal boundary layers, we can explain reasonably well the differences between the fully compressible and anelastic models, in terms of both the heat transfer and viscous dissipation dependence on compressibility. This could be mostly an effect of density variations on thermal diffusivity. Based on the different forms of entropy balance between exact and anelastic models, we find that a necessary condition for convergence of the anelastic results to the exact solutions is that the product $\unicode[STIX]{x1D716}q$ must be small compared to unity, where $\unicode[STIX]{x1D716}$ is the ratio of the superadiabatic temperature difference to the adiabatic difference, and $q$ is the ratio of the superadiabatic heat flux to the heat flux conducted along the adiabat. The same condition seems also to be associated with a convergence of the computed heat fluxes. Concerning the anelastic liquid approximation, we confirm previous estimates by Anufriev et al. (Phys. Earth Planet. Inter., vol. 152, 2005, pp. 163–190) and find that its results become generally close to those of the fully compressible model when $\unicode[STIX]{x1D6FC}T{\mathcal{D}}$ is small compared to unity, where $\unicode[STIX]{x1D6FC}$ is the isobaric thermal expansion coefficient, $T$ is the temperature (here $\unicode[STIX]{x1D6FC}T=1$ for an ideal gas) and ${\mathcal{D}}$ is the dissipation number.